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of the hadronic component of the photon is also studied. We examine the 
possibility of distinguishing between different parametrizations of the photon 
parton densities using charm production data, and the possibility of extracting 
information about the small- x behaviour of the gluon density of the proton. 
We also give a prediction for the transverse momentum and pseudorapidity 
distributions for bottom production at HERA. 
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1. Introduction 

Experimental results on charm photoproduction at HERA have recently become 
available J1J. It is likely that more detailed results on the differential distributions 
will appear in the near future. With respect to previous photoproduction experiments 
(ref. |2J) HERA offers the new opportunity of a higher energy regime. A theoretical 
study of the total photoproduction cross section has already been given in ref. [[J . In 
the present work we extend the analysis of ref. |§ (to which we refer the reader for a 
general introduction and for notation) by considering single inclusive distributions for 
heavy flavour production at HERA. Our analysis is based on the next-to-leading order 
calculation of heavy-quark photoproduction and hadroproduction cross sections per- 
formed in refs. as implemented in a computer program developed in refs. |||7[]. 
We will consider cross sections at fixed photon energies. When necessary, a discus- 
sion of the energy dependence of the distributions will be given. We will also consider 
electroproduction for small photon virtuality, in the Weizsacker- Williams approxima- 
tion, which is appropriate for the bulk of the photoproduction cross section. We will 
not consider the production of heavy flavours in deeply inelastic events, i.e. events 
in which high photon virtuality is required (see refs. ||). The paper is organized 
as follows. In section 2 we study the point-like contribution of the transverse mo- 
mentum and pseudorapidity distributions in photoproduction. Section 3 is devoted 
to the study of the hadronic component, considering the possibility of separating it 
from the point-like component with appropriate cuts. In sections 4 and 5 we consider 
electroproduction of charm and bottom, in the Weizsacker- Williams approximation. 
In section 6 we compare our fixed-order results with those of the Monte Carlo HER- 
WIG [[|, and in section 7 we give our conclusions. Some technical details on the 
factorization schemes, the scale dependence, and the modified Weizsacker- Williams 
approximation used here are given in the Appendix. 



2. point-like component 

Unless specifically stated, we will use in the following the set of parton densities 
MRSA, ref. ||10||, with A 5 = 151 MeV. This set of distribution functions has been 



recently updated (MRSG, ref. to include new HERA deep inelastic scattering 

data, allowing a different small- x behaviour of the gluon and sea quark densities. We 
have checked that the shapes of the single-inclusive distributions we are considering 
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are not significantly different when the new parametrization is used. The default 
value of the charm quark mass will be m c — 1.5 GeV. The renormalization scale will 
be taken as fi R = fiQ, and the factorization scale for the proton and for the photon 
will be taken = 2/xo, where /xo = \Jp% + m c- We begin by showing the 

transverse momentum distribution for the charm quark in fig. [I], for E 1 = 25 GeV 
(the proton energy E p will be fixed in the following to be 820 GeV). We also show 




Figure 1: Charm transverse momentum distribution in photon-proton colli- 
sions, with and without a pseudorapidity cut. The effect of applying a Peterson 
fragmentation function to the final-state quark is also shown. 



the effect of applying the Peterson fragmentation function 

D(x) = j (2.1) 

V ; x{l-l/x-e/{l-x)f V ; 

to the final state quark. We use the value e = 0.06, which is the central value 
quoted in ref. for charm quarks. As can be expected this softens considerably 
the p T spectrum. As discussed in ref. [|L4[ , in the case of fixed target photoproduction 
experiments, the inclusion of Peterson fragmentation correctly reproduces the shape 
of the measured p T distribution. The effect of a pseudorapidity cut similar to the one 
applied by the ZEUS collaboration |I[ is also shown. At low to moderate p T this cut 
considerably lowers the cross section, while at higher p T it has a negligible effect. In 
fig. El we compare the p T distributions at different photon energies. We also illustrate 
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the effect of applying an intrinsic transverse momentum kx to the incoming parton 
(see ref. We find that even with the very large value of (k T ) = 2 GeV 2 the 

effect is small. Therefore, in the following, we will neglect the effect of an intrinsic 
transverse momentum of the incoming partons. 
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Figure 2: Charm transverse momentum distribution in photon-proton colli- 
sions, for two values of the photon energy. The effect of applying an intrinsic 
transverse momentum to the incoming parton is also shown. 




Pseudorapidity distributions for charm are shown in fig. ||| for E 1 = 25 GeV. We 
observe that the point-like contribution to the cross section at this energy strongly 
favours large negative pseudorapidities. A cut on the transverse momentum of the 
produced quark tends to move the distribution towards the central region. The dotted 
curve shows the effect of fragmentation in the presence of a transverse momentum 
cut. The fragmentation has little effect on the pseudorapidity of the quark, but 
degrades its transverse momentum, so that the p T cut has a stronger effect. The 
effect of fragmentation without transverse momentum cut is not shown in the figure, 
since it is not well defined (see ref. ]14)). At lower energies (fig. |j) the pseudorapidity 
distributions become more central. 

We now turn to the sensitivity of our distributions to the various parameters that 
enter the computation. We have studied the dependence of our distributions upon 
the charm quark mass, by varying it between 1.2 and 1.8 GeV, on the renormalization 
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7P -> c+X, pointlike component 

E 7 =25 GeV, E p =820 GeV 

Solid: QCD NLO 

Dashed: QCD NLO, p T >2 GeV 

Dotted: QCD NLO + fragm., p T >2 GeV 




Figure 3: Charm pseudorapidity distribution in photon-proton collisions, 
with and without a transverse momentum cut. The effect of fragmentation is 
also shown. 
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Figure 4: As in fig. for E 1 = 3 GeV. 
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scale, which was varied by a factor of 2 below and above its default value, and upon 



the distribution functions, by considering the MRSD-' JT5| and the CTEQ2MF |IB| 
sets. The proton parton densities were chosen in order to span the allowed range for 
the small-x behaviour of the gluon density, MRSD-' being the most singular, and 
CTEQ2MF the least singular one (see ref. || for a discussion of this point). Studying 
the variation of our results with respect to A is a difficult problem, since the values 
of A extracted from fits to deep inelastic scattering data are not in good agreement 
with the LEP value. The CTEQ group has provided a parton density set in which 
the value of A was pushed as high as possible (see also ref. fl7j). They use a value 
of A 5 = 220 MeV, which is still below the LEP value, but adequate for a study of 
the sensitivity of the distributions. The set of parton densities obtained in this fit 
(CTEQ2ML) was used in association with this value of A. 

We found that the only cases in which there is a noticeable shape variation is in the 
p T distributions when we vary the mass and the parton densities. The relevant plots 
are given in figs. ^] and From the figures we can see a comparable variation of the 
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Figure 5: Mass dependence of the p T distribution of charm. 



p T spectrum in the low momentum region in the two cases. The mass dependence is 
easily understood. We expect nearly no mass dependence at high transverse momenta, 
while in the massless limit the cross section diverges at small momenta. Therefore, 
the smaller the mass, the higher the cross section at low p T . As far as the distribution 
function dependence is concerned, the small- a; uncertainty in the densities plays here 
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Figure 6: Sensitivity of the p T distribution of charm to the parton densities. 



a major role. The more singular the small- re behaviour, the higher the cross section 
at small pr- 

As a last point, we observe that, when measuring the total cross section, if a 
transverse momentum cut on the p T spectrum is being applied, the extrapolated 
total cross section will depend upon the assumptions we make about the value of the 
mass and the small-x behaviour of the distribution functions. This has to be properly- 
taken into account when inferring properties of the parton density at small x from 
the energy behaviour of the total charm cross section. 



3. Hadronic component 

The point-like contribution of the charm cross sections will be contaminated by 
the hadronic component. Depending upon the chosen distribution functions for the 
photon, the contribution of the hadronic component to the total cross section may 
dominate over the point-like one. Whatever the choice of the distribution functions, 
however, the two components differ remarkably in the pseudorapidity distribution. 
The point-like component, as we have seen, favours negative pseudorapidities, while 
the hadronic component favours the positive direction. The separation of the two 
components will therefore require an analysis of a large pseudorapidity range. In the 
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following, we will show a study of the full (i.e. point-like plus hadronic) pseudorapidity 
distribution for different choices of photon and proton distribution functions. We will 
see that the shape and magnitude of the rj distribution will allow us to distinguish 
between the most extreme parametrizations of photon distribution functions. In 
fig. we show the r] distribution, obtained with the LAC1 ]18[ set of photon parton 
densities, and three different sets of proton parton densities; both the point-like and 
the full results are presented. A similar plot for the GRV |0| set for the photon is 




Figure 7: Pseudorapidity distribution of charm quarks, obtained with the 
LAC1 set of photon parton densities, and three different sets of proton parton 
densities. The point-like component is also shown. 



given in fig. ||. We have chosen the LAC1 and GRV sets because, for our purposes, 
they represent the two extreme possibilities, as discussed in ref. [§]. In figs. [?] and 
[8] we have imposed a realistic transverse momentum cut, which in general will make 
the pseudorapidity distribution even narrower. In spite of the cut, the separation of 
the two contributions is quite apparent. 

The study of the pseudorapidity distribution could help in distinguishing among 
different proton parton density sets, especially if the large negative pseudorapidity 
region could be explored. On the other hand, the large difference in shape induced 
by the two photon sets considered has measurable effects even in the central region. 

The full transverse momentum distribution is presented in fig. |[ The point-like 
and hadronic contributions are also separately shown. From the figure, we can see 



-8- 




-4 -2 2 4 

n 

Figure 8: Pseudorapidity distribution of charm quarks, obtained with the 
GRV set of photon parton densities, and three different sets of proton parton 
densities. The point-like component is also shown. 

that the hadronic contribution is dominant in the low-p T region. At p T ~ 2 GeV, it 
amounts to 50% of the full cross section. Therefore, even without investigating the 
very low-p T region, the softer behaviour of the full cross section with respect to the 
point-like one should be visible. Obviously, this behaviour is strongly influenced by 
the photon parton densities used. If the GRV set were used, instead of the LAC1 
set, the full distribution would differ only sligthly from the point-like one already 
at moderate p T . We can conclude that the charm transverse momentum spectrum 
at large values of the photon energy could help in further constraining the parton 
densities in the photon. From fig. [7], we also learn that asymmetric pseudorapidity 
cuts may enhance or suppress the effect of the hadronic component in the p T spectrum, 
thus providing another handle for its study. 

4. Charm electroproduction 

In this section we discuss single-inclusive distributions for charm production in 
electron-proton collisions in the Weizsacker- Williams approximation. We take the 
electron and proton energies to be 27.5 GeV and 820 GeV, respectively. Since the 
Weizsacker- Williams function peaks in the small- X region, the bulk of the contribution 
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Figure 9: Transverse momentum distribution of charm quarks. The point- 
like and hadronic contributions are separately shown, together with the sum of 
the two. 



to the cross section will be due to photons of relatively low energy. Therefore, the 
role of the hadronic component is less important with respect to the monochromatic 
photon case. By applying a small-p T cut, the contribution of the hadronic compo- 
nent is even more suppressed, especially for soft photon parton densities. Therefore, 
the QCD predictions at next-to-leading order are less affected by the uncertainty 
originating from the photon distribution functions. 

Because of the softness of the incoming photons, the pseudorapidity distribution 
of the point-like component is concentrated in the central region. This can be seen 
in fig. ITDL where we show the pseudorapidity distribution for charm quarks, supple- 
mented with Peterson fragmentation and a transverse momentum cut. The hadronic 
component is computed using our two extreme sets of photon parton densities. It 
is apparent that its effect is small in the central pseudorapidity region explored at 



HERA. The curves in fig. |10| are obtained with m c = 1.5 GeV, fx R = /io, /i F = /i 7 = 2/i 
and the MRSA set for proton parton densities. While the total cross section is quite 
sensitive to the value of these parameters, as shown in ref. |3[], the shape of the single- 
inclusive distributions is rather stable. The pseudorapidity distribution shape is only 
mildly dependent upon the choice of proton parton densities. This is due to the fact 
that the various parametrizations differ in the small-x region, which is not deeply 
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Figure 10: Pseudorapidity distribution for charm electroproduction. The 
proton parton density set MRSA is used. 



probed in the electroproduction process. This is shown, for the point-like compo- 
nent, in fig. |TT[ The proton parton density dependence of the hadronic component is 
completely negligible, due to the softness of the partons in the electron. 




Figure 11: Proton parton density dependence of the pseudorapidity distribu- 
tion for charm electroproduction; only the point-like component is shown. 
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We now consider the transverse momentum distribution. In fig. [12| we show the 
point-like contribution to this differential cross section, with and without Peterson 
fragmentation. The softening of the distribution due to fragmentation should be 
observable at HERA. The effect of applying a pseudorapidity cut is also shown: it 
affects the shape of the distribution less dramatically than in the case discussed 
in section 2. We have checked that the hadronic component contribution to this 
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Figure 12: Transverse momentum distribution for charm electroproduction 
(point-like contribution). 



distribution is remarkably softer than the point-like one. In practice, its effect is less 
than 10% for p T > 2 GeV when fragmentation is included and the pseudorapidity cut 
| T] | < 1.5 is applied. We have also computed the same distribution for different values 
of the renormalization and factorization scales, and we found that the corresponding 
shape variations are small. 

All the distributions presented so far were also evaluated in the case when an 
antitag condition on the outgoing electron is applied, that is to say when the elec- 
tron scattering angle is below a given value, which we chose to be 5 mrad. For this 
computation, the improved form of the Weizsacker- Williams distribution presented 
in ref. [^(J is appropriate (see the appendix). Differences in the shapes of the distri- 
butions are found to be small. As an example, we show in fig. ^ the pseudorapidity 
distribution with and without the antitag condition described above. The antitag 
condition has the effect of decreasing the total cross section, as described in ref. 0], 



of slightly enhancing the contribution of positive pseudorapidities, and of mildly soft- 
ening the transverse momentum spectrum. 




Figure 13: Pseudorapidity distribution for charm electroproduction, with and 
without an antitag condition on the scattered electron. 



We can conclude that charm electroproduction distributions are rather insensi- 
tive to the choice of the parameters entering the perturbative calculation. Therefore, 
single-inclusive charm electroproduction will be of little help in constraining the par- 
ton distribution functions of both the proton and the photon. On the other hand, 
the comparison between data and theoretical predictions will be useful for the study 
of the production mechanism. 



5. Bottom electroproduction 

Due to the higher value of the quark mass, perturbative QCD predictions for bot- 
tom production are more reliable than those for charm. Furthermore, we are able to 
study the scale dependence in an exhaustive manner. In fact, the factorization scale, 
when varied by a factor of two below and above its default value, is always higher 



than the mininum allowed by the distribution function parametrizations. In fig. [14 
we show the point-like component of the pseudorapidity distribution for bottom elec- 
troproduction. The solid curve represents our central prediction, which corresponds 
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to choosing m& = 4.75 GeV, \l r — \i F — fi 7 — fio, where fio = yp^ + ml. The QCD 
result is supplemented with Peterson fragmentation with e = 0.006, the central value 
of ref. [13| for bottom quarks. In the case of bottom the effect of fragmentation is 
less important than in the case of charm. The band between the two dashed curves 
has been obtained by varying the bottom quark mass between 4.5 and 5 GeV and the 
renormalization and factorization scales between /zo/2 and 2/io- The proton parton 
density set chosen is MRSA. We also show (dotted curve) the result obtained with a 
parton density set fitted at a higher value of A QCD , namely the set CTEQ2ML with 
A 5 = 220 MeV, and the central values of mass and scales. It is difficult to foresee 




Figure 14: Pseudorapidity distribution for bottom electroproduction (point- 
like component only), with Peterson fragmentation and a transverse momen- 
tum cut. 



what a realistic transverse momentum cut could be for bottom production, and we 
have therefore kept the same cut p T > 2 GeV we used in the charm case. 

Figure |T4| does not describe the uncertainties on the theoretical prediction com- 
pletely, because the CTEQ2ML set was used keeping the scales and the mass at their 
default value, and because the hadronic contribution was not included. The full un- 



certainty on the 7] distribution is presented in fig. [15], in which mass and scales are 
varied also when using the CTEQ2ML set, and the possibility is considered that the 
effect of the hadronic component is as large as the LAC1 set implies, or is completely 
absent. The dramatic effect of the inclusion of the hadronic component is apparent 
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Figure 15: Full uncertainty on the pseudorapidity distribution for bottom 
electroproduction with Peterson fragmentation and a transverse momentum 
cut. 



in the positive pseudorapidity region. This is due to the very soft behaviour of the 
LAC1 partons and to the stronger sensitivity of the hadronic component to the choice 
of mass and scales with respect to the one of the point-like component. 



We now turn to the transverse momentum distribution, shown in fig. |T(J, which 
is analogous to the corresponding figure for the 77 distribution, fig. [14]. As expected, 
in the large-p T region the sensitivity to mass and scales choice is strongly reduced. 
Therefore, the perturbative prediction becomes more reliable there. 



The full uncertainty on the p T distribution is given in fig. [T7|, which is analogous 
to fig. [TR From the figure, it is quite clear that even with the LAC1 set the hadronic 
component affects the prediction only marginally. This fact is a consequence of the 



applied pseudorapidity cut, as can be inferred from fig. [L5|. We can therefore regard 
fig. [H] as a reliable prediction of QCD for the p T spectrum of b hadrons at HERA. 



6. Comparison with Monte Carlo results 



In order to assess the influence of higher-order and non-perturbative contribu- 
tions on the shape of single-inclusive distributions, we have compared our fixed-order 
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Figure 16: Transverse momentum distribution for bottom electroproduction 
(point-like component only) with Peterson fragmentation and a pseudorapidity 
cut. 
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Figure 17: Full uncertainty on the transverse momentum distribution for 
bottom electroproduction with Peterson fragmentation and a transverse mo- 
mentum cut. 



-16- 



calculations with the results obtained with the parton-shower Monte Carlo program 
HERWIG 0. When using HERWIG we have always excluded the flavour excitation 
processes. This is because comparable processes, such as gluon splitting, are diffi- 
cult to include. We have therefore preferred to use the Monte Carlo at a consistent 
level of accuracy. As discussed in section 2, the effect of an intrinsic transverse mo- 
mentum of the incoming partons is negligible in photoproduction at HERA energies. 
This leads us to expect that higher-order and non-perturbative corrections are not as 
important here as in the case of hadropro duct ion [14]. On the other hand, it is inter- 
esting to investigate how Peterson fragmentation compares with the model of cluster 
hadronization implemented in HERWIG. We limit our discussion to the point-like 
component of the photoproduction cross section for fixed photon energy. 

In fig. pl| we show the transverse momentum distribution computed in fixed-order 
perturbation theory and with HERWIG. All HERWIG curves have been normalized 
to give the same total cross section of the corresponding next-to-leading order QCD 
predictions. We observe that before the hadronization process is switched on, the 
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Figure 18: Transverse momentum distribution for charm photoproduction 
(point-like component) as given by next-to-leading order QCD and by the 
Monte Carlo HERWIG. 



QCD prediction and the HERWIG result agree almost perfectly in all the presented 
p T range (upper curves). We then consider the effect of hadronization in the central 
pseudorapidity region (\r}\ < 1.5) by comparing the perturbative prediction supple- 
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mented with Peterson fragmentation and the HERWIG results for charmed hadrons 
or D mesons only (lower curves). In this case the two predictions display a different 
behaviour only in the low-p T region. This is also the only part of the p T spectrum 
that is significantly affected by the pseudorapidity cut. The results in fig. |Tj] allow us 
to conclude that Peterson fragmentation gives a description of the hadronization pro- 
cess, which is consistent with the one suggested by cluster models only at sufficiently 
large transverse momenta. 



A similar study for the pseudorapidity distribution is presented in fig. 19. In this 
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Figure 19: Pseudorapidity distribution for charm photoproduction (point-like 
component) as given by next-to-leading order QCD and by the Monte Carlo 
HERWIG. 



case, some difference is observed between the two distributions before hadronization, 
the fixed-order QCD prediction being slightly broader. However, they are peaked 
around the same value of rj, which is rather far from the central region, where the 
experimental acceptance is larger. Applying a small transverse momentum cut, the 
use of Peterson fragmentation becomes possible also for the pseudorapidity distribu- 
tion. We are therefore able to compare the QCD next-to-leading order distribution, 
convoluted with the Peterson fragmentation function, and the HERWIG curves for 
charmed hadrons. We find a complete agreement between the two. 

We have performed the same kind of analysis on the hadronic component, and 
we found that the qualitative behaviour is completely analogous to what has been 
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presented for the point-like component. 

7. Conclusions 

We have studied single-inclusive cross sections for charm and bottom production at 
HERA. We considered charm production in monochromatic photon-proton collisions, 
for energies spanning the range accessible at HERA. We found that the shape of the p T 
and pseudorapidity distributions are quite stable with respect to the variation of the 
input parameters entering the calculation. On the other hand, the total cross section 
value (see ref. J| for a discussion on this point) is strongly affected by the choice 
of parameters. We found that from the study of the pseudorapidity distributions 
in the central rj region (\rj\ < 1.5) it is difficult to distinguish between the various 
proton parton densities with different small- x behaviour. On the other hand, the 
striking difference between the predictions obtained by using the LAC1 and GRV set 
for photon distribution functions is clearly visible. This fact is almost completely due 
to the very soft behaviour of the LAC1 parton densities. 

The transverse momentum distribution for large p T is reliably predicted by QCD. 
Measurements in this region should constitute a good test of the production mech- 
anism. At moderate p T , the point-like result is strongly influenced by the hadronic 
contribution when the LAC1 set for photon parton densities is used. Experimental 
investigations in this region of the spectrum could therefore help to distinguish among 
the various parton distribution functions in the photon. 

Charm electroproduction is a less clean test of the production mechanism with 
respect to monochromatic photon production, because of the convolution with the 
Weizsacker- Williams function. On the other hand, the softness of the latter implies 
a strongly reduced dependence of the total cross section upon the phenomenological 
parameters entering the calculation. Furthermore, when a small transverse momen- 
tum cut is applied, the hadronic contribution is more suppressed than the point-like 
one, whatever the photon densitites used. Therefore, electroproduction is less sen- 
sitive to the contamination of the hadronic component than photoproduction with 
a monochromatic photon, and can be used to further investigate the point-like pro- 
duction mechanism. Unfortunately, a sizeable dependence upon charm mass and 
renormalization scale remains in the total cross section prediction. 

The case in which an antitag condition is applied on the outgoing electron was 
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also studied. Aside from the effect on the total cross section value (which has already 
been discussed in ref. ||), the angular cut slightly enhances the importance of the 
central pseudorapidity region and mildly softens the transverse momentum spectrum. 

We also considered the production of bottom quarks in ep collisions. Due to the 
higher value of the quark mass, perturbative QCD predictions for bottom production 
are more reliable than those for charm. It turns out that the point-like component 
is quite stable in shape. The hadronic component has a dramatic effect on the pseu- 
dorapidity distribution, although less important than in the case of charm. Its effect 
in the central rj region is however marginal. We are therefore able to give a reliable 
prediction of the p T spectrum for b hadrons in the central 77 region. 

Finally, we compared our results with those obtained using the parton shower 
Monte Carlo HERWIG. We limited our discussion to the point-like component for 
monochromatic photon-proton collisions. As far as open charm quark is concerned, 
we found perfect agreement for the p T distribution, in the whole range considered. 
The next-to- leading order QCD prediction for pseudorapidity distribution is instead 
slightly broader than the one given by HERWIG in the region around the peak. This 
difference is however immaterial in the central pseudorapidity region, which is ex- 
plored by the experiments. We also point out that, without any small transverse 
momentum cut, differential cross sections in "longitudinal" quantities such as rj are 
indeed not characterized by a true hard scale. In this case, differences between pertur- 
bative QCD and parton shower approach are expected. Coming to charmed hadron 
production, we compared HERWIG prediction with QCD calculation supplemented 
with the Peterson fragmentation function. A small-p T cut was applied. The agree- 
ment is surprisingly good in the whole r\ range, while the two predictions for transverse 
momentum distribution differ in shape only in the small-p T region. 
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Subtraction schemes in photoproduction 

A differential photoproduction cross section can be written as 

d*lT\P 7 , P H ) = da££(P 7 , P H ) + da£g(P 7 , P H ) , (.1) 

where the quantities in the right-hand side of this equation are the so-called point-like 
(or direct) and hadronic (or resolved) photon cross sections. In QCD 



^^point(-^7> P h) — / dxfj \x, /j,]r)dcr<yj(xP H: q; s (/x J j), /j, r , p, F , // 7 ) (.2) 

rf(7 hIdr( P 7; p h) = Yl j dxd yfi l] (xi lh)fj H \v> ^F)dcr tj (xP y , yP H , a a (p, R ),fi B , /x 7 ) . 

(.3) 

In the case of heavy quark production at next-to-leading order, the short distance 
partonic cross sections are given by 

+a em a^(/i H )(i(T 7 j ) (pi, p 2 , A*h, A*f, A* 7 ) (-4) 
da ij {p 1 ,p 2 ,a s {fi R ),fi R ,fi F ,fi 1 ) = a 2 s (/i B )da>f(p 1 ,p 2 ) 

+al(fi R )dalj } (pi,p 2 , A* fl , //f, /x 7 ) • (-5) 

In dimensional regularization, we have 
dcr^(p 1 ,p 2 ) = da$(pi,p2, ~) + — Y,J dx (zPk 7 {x) ~ H kl (x)j daf] {xp x ,p 2 ) 



+7T~ J dx (= p kj( x ) ~ K kj\ x )) da^kiPu xp*)) (-6) 

k 

d^f(pi,P2) = daff(p 1 ,p 2 , -) + — J2f dx (= p ki( x ) ~ K kh x )) dcr i ^{xp 1 ,p 2 ) 
+ / dx(±P kj (x)-Klf(x)\da%\p u xp 2 ), (.7) 

k 

where da^zj (dcr^)) and da^ 1 (daff) are the full, <i-dimensional regulated partonic 
cross sections at the next-to-leading (leading) order for the point-like and hadronic 
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contributions respectively. The 1/e singularities in daff and da^- are appropriately 
subtracted on the right-hand side of eqs. (0), ([7|). 

The Pij are the usual Altarelli-Parisi splitting function; P J7 is the kernel for the 
splitting process 7 — > j +j, which is equal to P gg up to a colour factor. The functions 
Kj?\ K^J and are completely arbitrary, in that they define an extra finite part 
of the subtraction; different choices correspond to different subtraction schemes. In 
eqs. ([TBj) and ([7]) the MS scheme is equivalent to H = K = 0. For greater generality 
we have admitted the possibility to have different subtraction schemes on photon and 
hadron legs. A change in the subtraction scheme implies that parton distribution 
functions are modified as follows 

// = /i + ^* + ^E^®/i- (-8) 
The term H is present only in the photon case, and it is a direct consequence of the 



inhomogeneous term in the modified Altarelli-Parisi equations. In ref. |22] a factor- 
ization scheme (DIS 7 ) for the photon densities is introduced, which uses K — and 
H 7^ 0. As can be seen for eqs. ([6|) and (0), the DIS 7 scheme is therefore equivalent 
to MS as far as the hadronic component is concerned. It does however modify the 
point-like component through the second term in the right-hand side of eq. (|1|), which 
is also responsible for a subtraction on the photon leg. This clearly shows that the 
point-like and the hadronic component are closely related, and that they should not 
be considered separately. In principle, when considering the variation with respect to 
/x 7 , cancellations occur between the hadronic and the point-like components, resulting 
in renormalization group invariance of the physical cross section da^ lH ^ up to terms 
of order a em a|. 

Photoproduction cross sections are usually linked to electroproduction ones by a 
convolution with the Weizsacker- Williams function |23|| : 

(.9) 



d^ eH \P e ,P H ) = Jdxf(x,fi ww )da^ H \xP e ,P H ). 



At HERA, the HI and ZEUS experiments can directly investigate the virtuality of the 
exchanged photon, by tagging the emitted electron and retaining only those events in 
which the electron scattering angle 9 satisfies the condition 9 < 9 C , with 9 C typically 
of the order of few mrad. In ref. |2(J a functional form especially suited for this 
experimental set-up was proposed, namely 



f{x,E e 



0., 



2tt 



2(1 



x) 



m 2 e x 



Em 



x) 2 9% + mix 2 



x 
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l + (l-x) 2 El(l-xYel + mlx* \ 

h ~ ° g -~,2~.2 | ' V ' 



x mtx* 



where E e is the incoming electron energy in the laboratory frame. When no angular 
cut is applied, eq. ([9|) is too restrictive; we relax it by writing 



du^ \P e , Ph) — J ' dxj ' {x^ ww )dcF^^ t {xP e) P H ) + J dxf(x, fx ww )do'l^ a J,( y xP e , P H 



.11^ 



with fi ww ^ n' ww (for a discussion of the \i ww and \i! ww values, see ref. and 
references therein). Strictly speaking, the fact that \i ww ^ jj,' ww gives up the renor- 
malization group invariance of the physical cross section in eq. (TT), in that no 



cancellation occurs between the point-like and the hadronic component when varia- 
tion with respect to /x 7 is considered. In practice, we have verified that this effect is 
numerically negligible, the point-like component being almost invariant with respect 
to /i 7 variations. 
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